Developed Smectics: When Exact Solutions Agree 
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In the limit where the bending modulus vanishes, we construct layer configurations with arbitrary 
dislocation textures by exploiting a connection between uniformly-spaced layers in two dimensions 
and developable surfaces in three dimensions. We then show how these focal textures can be used 
to construct layer configurations with finite bending modulus. 



When subject to frustrating boundary conditions or 
extreme strains, liquid crystals, superfluids, and mag- 
nets will locally rise into their higher-symmetry phases 
resulting in point, line, and planar defects pQ. Energetic 
considerations determine the dimensionality of these de- 
fects; in some systems, rigorous results demonstrate that 
the energy minimizers will have point or line defects 
Smectic liquid crystals represent a special challenge as 
they are described by an essentially nonlinear elasticity 
theory [TJ HHB] that gives rise to anomalous elasticity [7J , 
dynamics [8 , 9J , and qualitatively modified ground states 
[Tffl ITT] . These nonlinearities are generic features of elas- 
tic systems with free surfaces [T3] and, thus, smectics are 
ideal systems for understanding elastic geometric nonlin- 
earities in general. Previously, we have studied smectic 
liquid crystals in the limit where the bending energy is 
neglected so that the layer spacing is strictly constant 
[T51 ITT] . Here we extend some of these techniques by 
employing a connection between developable surfaces in 
three dimensions and uniformly-spaced layers in two di- 
mensions, allowing us to find layer configurations for any 
specified dislocation texture. We compare these solu- 
tions with exact solutions to the nonlinear elasticity [10] 
equations which only allow superposition of dislocations 
along a single line [TBI, ITrj] . Not only do the two solu- 
tion methods agree, but the geometric construction ex- 
plains the fundamental asymmetry of the smectic strain 
field around a dislocation, first predicted by Brener and 
Marchenko [10], and sheds light on the simple, topo- 
logically based, Bogomornyi-Prasad-Sommerfield (BPS) 
bound [16] [IT]. We exploit this understanding to con- 
struct textures for dislocations separated by a finite num- 
ber of layers with finite bending rigidity. 

The order in a smectic is characterized by the phase 
field 0(x) appearing in the density modulation bp oc 
cos[27T0(x)/a], where a is the natural layer spacing. In 
terms of <f> the free energy is the sum of compression and 
bending contributions 
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where B is the compression modulus, A 
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is the penetration length and Ki is the bending mod- 
ulus. In smectics A, the normal to the smectic lay- 
ers is the nematic director n = V0/|V0|. Geometri- 
cal and topological insight is gained by considering the 
surface [x, y, 4>(x, y)\ £ R 3 with surface normal N = 



[-d X <t>,-dy<f>, 1] /V 1 + |W| 2 PI- Here > WC? ShaU foCUS 

our attention on the limit A < a, or K\ — ¥ 0, where 
bending becomes unimportant compared to compression. 
Indeed, when A = the free-energy is strictly minimized 
when |V</>| = 1; differentiating (V<^) 2 = 1, we have: 
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which requires the Gaussian curvature, K oc d 2 (f>dy(f> — 
(d x d y <fi) 2 = 0. It follows from GauB's Theorem Egregium 
that our surface must be isometric to the plane, so it can 
be built out of sections of planes, cones, cylinders, and 
tangent-developable surfaces. The constant angle condi- 
tion further restricts to planes, cones, and the develop- 
ment of cylindrical helices [TO] . 

It is amusing that the latter can be used to gener- 
ate uniformly-spaced involutes of curves j^OJ [2JJ ; though 
known to the ancients [22] , we will briefly review the 
connection between level sets of constant-angle, devel- 
opable surfaces and involutes. Consider a curve R(cr) = 
[x(a),y(a), z(a)] in K 3 , parameterized by its arclength 
<7, with Frenet-Serret frame [t, v,0\ = [R.,t/n,t x is], 
curvature n(a) > and torsion r. The tangent de- 
velopable surface is defined in terms of the curve and 
its family of tangents: X(<7i,<72) = R-(ci) — ctit(o\) 
for a 2 > 0. Note that the unit normal to the sur- 
face N(cti,ct 2 ) = <9iX x a 2 X/|5iX x <9 2 X| = f3(o x ), the 
curve's binormal at <j\. It follows that N only depends 
on <j\ and so the Gaufi curvature vanishes. If the angle 
between N and z is constant then so is the angle be- 
tween j3 and z. Differentiating with repect to a, we have 
= z • = —tz ■ v so v lies in the xy-plane. Define the 
surface curve ~f(s) — X(cri+s, (72+s) with tangent ~f(s) = 
(cr 2 + s)k((Ti + s)v(ai + s). 7 lies in a plane of constant 
z = c and y(s) sweeps out an involute starting at s = on 
the planar curve R_l(s) = [x{a\ + s), y(<J\ + s), c]. Apart 
from concentric circles, any set of uniformly-spaced invo- 
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FIG. 1: (Color online) We construct the two-dimensional lay- 
ers by taking level sets of a piecewise developable (Gaussian 
curvature K = 0) surface which makes a constant angle with 
the z direction. From back left to close right, the surface 
is made of two intersecting planes which end and attach to 
pieces of cone which necessarily intersect another set of par- 
allel planes on a parabola. 



Taking level sets of the surface produces an uniformly- 
spaced smectic texture for a dislocation. Aside from the 
point defect corresponding to the cone's vertex, there 
is a focal set consisting of the two parabolic segments 
x 2 = b\y\ + 6 2 /4. Recall that in the linear theory the elas- 
tic response is concentrated in two full parabolic regions 
above and below the defect [23]. The present construc- 
tion only generates compression strain on the "right" side 
of the defect. 

In the presence of a defect, BPS minimizers of ([T]) and 
related free energies were found [TO] [HI HH] and, for small 
^/y [H]j the displacement for a single defect at (x,y) = 
(0, 0) is u(x, y) = y- cj>(x, y) is 



u(x,y) = 2Asgn(y) In 
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where E(x) = (tt) 1 / 2 J* dtexp(— t 2 ) is the error func- 
tion. The associated compression strain e for y > scales 
as 
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lutes will generate an evolute curve which constitutes a 
singularity or edge of the surface and where the bending 
of the involutes diverges. Since this will generate a two- 
dimensional region without smectic order, we will not 
consider such cases, although surfaces like this are liable 
to play a role in sample cells with large inclusions. Here 
we are interested in defects that can be reduced to points 
and lines and so we only consider constant angle cones 
and planes. For convenience we set the constant angle to 
be 7r/4. 

Smectics enjoy two types of point defects, disclinations 
and dislocations. In the language of surfaces the discli- 
nations are critical or singular points on the graph of 
4>. Dislocations can be constructed by choosing </> = 
x + (6/27r) arg (x + iy) to be a tilted helicoid [T5], re- 
sulting in a two-dimensional smectic with bending and 
compression deformations. However, we can also build a 
dislocation with vanishing compression with lines across 
which the director jumps discontinously, thus being vis- 
ible under light microscopy. To this end, consider the 
construction of an edge dislocation shown in Fig. [1] Two 
planes meeting along a ridge are connected to two sim- 
ilar planes, that meet along a ridge at a lower height 
(6/2 lower where 6 6 aZ is the Burgers scalar), by a por- 
tion of a cone. The cone's apex coincides with the end- 
point of the upper ridge and the transition from plane 
to cone is Lipschitz C . However, the intersection with 
the lower pair of planes introduces a cusp, or focal curve, 
along which the normal changes discontinuously and the 
surface is only Lipschitz C°, as is the director field. 
This focal curve consists of part of a pair of parabolas. 



For large A/6 this reproduces the symmetric, linear strain 
field. However, as A/6 — > 0, we have d y u ~ 9(x)S(y — 
x 2 /b), half of a parabola on the side with fewer layers, 
and, as shown in Fig. [2] in agreement with the focal con- 
struction. Though the shape of the parabola is identical 
in the focal and BPS solutions, we note that there is a 
vertical offset of 6/4 between them. Because the BPS 
solution is based only on a step-function boundary con- 
dition at y = used to satisfy the topology at infinity, we 
do not expect the near-defect details to be reproduced, 
but for large x and y, the solutions agree as shown in 

Why should the strain be asymmetric [TD] ? Recall 
that the nonlinear compression strain e measures the de- 
viation of the wavenumber q = 2ir/d from qo = 2n/a, 
e oc (q — qo) 2 and so, away from the linear regime, com- 
pression d < a is more energetic than dilation d > a. It 
follows that in the equal-spacing limit, the texture will 
preferentially distort on the dilated side. The presence 
of a focal line in the BPS solution also is not a mystery. 
Differentiating the BPS equation 
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with respect to x yields the Burgers equation d y v — 
vd x v = Xd 2 v for v = d x u. As is well known, the in- 
viscid Burgers equation has straight characteristics and 
produces asymptotically parabolic shock curves as we 
have here [24 . In comparison, the focal construction 
arises from constructing characteristics of the geodesic 
condition (n • V) n = [13 . Expanding this equation to 
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FIG. 2: (Color online) Comparison of focal and BPS con- 
structions: the dashed lines show the layers from the fo- 
cal construction in Fig. [I] for the layers on one side of (a) 
one or (b) two dislocations. The solid curves are level sets 
of (j> = y — u(x,y), where it is determined by BPS evolu- 
tion, starting with the first layer of the focal construction, for 
A = 0.05 and A = 0.1. The background is shaded according to 
the compression energy of the asymptotic BPS solution, Eq. 
II]), for A = 0.05. The parabolic focal line is y = ±{x 2 /b-b/4); 
we use the vertical offset in the BPS solution [101116] . (c) The 
focal construction for two dislocations, built by attaching two 
single dislocations as in Fig. [I] Note that now some of the 
focal lines arise from the intersection of cones with cones and 
are pieces of hyperbolae, not parabolas. 



quadratic order in 6n s» n — y precisely yields Burgers 
equation in v — — x • <5n. 

When multiple defects lie along a line of constant y, 
the BPS method allows the superposition of defects via 
the Hopf-Cole transformation. We can superpose in the 
focal construction too: multiple edge dislocations can be 
constructed by repeating the procedure described for Fig. 
[T] For example, in Fig. [2] we show the construction for 
a pair, both located at the same value of y. Note that 
there are now new features: in addition to parabolic fo- 



cal curves, there are regions of the surface where cones 
intersect cones and, by definition, this happens along hy- 
perbolae. As we show in Fig. |3j it is also possible to con- 
struct arbitrary focal textures in which the dislocations 
no longer lie at the same value of y: when lines meet cir- 
cles they intersect on parabolae, when circles meet circles 
they intersect on a hyperbola. 

How does the BPS solution fare? Again we begin with 
the deformations for large \/y, where straightforward nu- 
merical analysis shows that hyperbolae are in the strain 
field 
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corresponding to a pair of dislocations |16j . Indeed, Fig. 
[2] shows remarkably good agreement between the focal 
construction and the BPS solution, including the details 
of the hyperbolae and the merging of the two focal curves. 

We also compute 'exact' solutions for the level sets 
4>{x,y) = y — u(x,y), shown as dark solid and dashed 
lines in Fig. [2j where the initial condition u(x, 0) is given 
by the phase field at y — in the focal construction. As 
expected from the asymptotic solution, the BPS evolu- 
tion respects the parabolic cusps in the focal construc- 
tion, deforming most to the left of the cusps but not 
on the right. This is to be expected; the deformation 
preferentially smooths out the higher curvature side and 
spreads the strain 'inside' the parabolic region in agree- 
ment with the predictions of linear elasticity. In BPS 
evolution, the quantities S± — exp{±u/(2A)} satisfy the 
extremal equations d y S± — ±\d%S± [THOU]. The evolu- 
tion has an inherent directionality: BPS evolution relaxes 
5+ to flat layers above the dislocation and SL below the 
dislocation. Therefore, a dislocation at y — requires 
the BPS evolution to change directionality on cither side 
of the line at y = 0. Similarly, it is possible to find the 
textures generated by multiple dislocations, as long as 
they lie along the y-axis. 

When defects sit at different values of y, we have to be 
more careful when A > 0. It is instructive to consider the 
difficulty in detail. First, consider the focal construction 
shown in Fig. [3] (dashed lines). In the vicinity of each dis- 
location, we expect the solutions at finite A to be approx- 
imated by BPS evolution. Above and below both dislo- 
cations, there is no difficulty constructing a valid BPS 
evolution since the BPS evolution directions agree. The 
layers between the two dislocations, however, must evolve 
upward on the left and downward on the right. We can 
reconcile this discrepancy by noting that the parabolic 
cusp between dislocations in the focal texture also forms 
a natural division between upward and downward evo- 
lution. As shown in Fig. ([3]), we evolve upward on the 
left of the parabolic cusp using the displacement for the 
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FIG. 3: (Color online) Two edge dislocations at arbitrary 
(xi,yi). The dashed curves are the focal layers, and the solid, 
gray curves are BPS evolution with A = 0.05. Parabolic (red 
and orange) and hyperbolic (green) cusps of the focal con- 
struction are also shown. The thick solid lines indicates the 
division between upward and downward evolving BPS solu- 
tions. 



lower dislocation as the initial condition. On the right, 
we evolve downward using the phase field for the upper 
dislocation as an initial condition. The result of evolv- 
ing upward and downward as indicated by the arrows 
in Fig. ([3]) is shown as solid layers. The layers arising 
from BPS evolution of opposite directionality meet nat- 
urally at the parabolic cusps without further adjustment 
because the deformation field is strongly asymmetric, in 
this case confined to the left of the parabolic cusp. Were 



this not to occur, we could, of course, impose continu- 
ity of the layers at the cusp by setting the displacement 
of the upward evolution equal to that of the downward 
evolution. The success of the focal method hinges on the 
asymmetry of the distortion field for small A. Once we 
have constructed the shape of the layer on either side of 
the two dislocations, we may continue the evolution out 
to infinity. Again, the BPS evolution preserves the un- 
derlying structure of the cusps of the focal textures and 
the regions of maximum strain (and layer deviation) oc- 
cur just to the left of the cusps. As long as the defects are 
further apart than A, this procedure should be reliable. 

In summary, we have developed a focal construction for 
multiple (and arbitrary) configurations of dislocations in 
a smectic. This construction uncovers a deep relation- 
ship between the BPS evolution of single and multiple 
dislocations and the focal construction. Using the nat- 
urally occurring cusps in the focal construction, we are 
able to develop BPS solutions for dislocations with layers 
between them that account for the geometric nonlinear- 
ities in the elastic strain. 
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